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(1.1) 



Abstract. In this paper, we prove the global existence and uniqueness of solution to d-dimensional 
(for d = 2, 3) incompressible inhomogeneous Navier-Stokes equations with initial density being 
bounded from above and below by some positive constants, and with initial velocity uo £ H"(R^) 
for s > in 2-D, or ito £ //^(R'') satisfying HitolLa [|Vuo||£2 being sufficiently small in 3-D. This in 
particular improves the most recent well-posedness result in [10] , which requires the initial velocity 
uo G H (R ) for the local well-posedness result, and a smallness condition on the fluctuation of 
the initial density for the global well-posedness result. 
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1. Introduction 

In this paper, we consider the global existence and uniqueness of the solution to the following 
d-dimensional (for d = 2,3) incompressible inhomogeneous Navier-Stokes equations with initial 
density in L°°(R'^) and having a positive lower bound: 

dtp + u-Vp = 0, (t, x) G R+ X R*^, 
p(dtu + u ■ Vu) — Au + Vp = 0, 
divn = 0, 

{p,u)\t=o = {po,uo), 

where p, u stand for the density and velocity of the fluid respectively, p is a scalar pressure function, 
and the viscosity coefficient is supposed to be 1. Such a system describes a fluid which is obtained 
by mixing two miscible fluids that are incompressible and that have different densities. It may also 
describe a fluid containing a melted substance. One may check [17] for the detailed derivation of 
this system. 

Given < po & L°°(R'^), and uq satisfying divuo = 0, ^/poUo G L^(R'^), Lions [17] (see also 
[5, 19] and the references therein for an overview of results on weak solutions of (1.1)) proved that 
(1.1) has a global weak solution so that 

^\\y/mu{t)\\h + ^ \\^u{r)\\l, dT < ^\\^uo\\l2. 

Moreover, for any a and f3, the Lebesgue measure 

G R"'; a < p{t, x) < /? | is independent of t. 

In dimension two and under the additional assumption that po is bounded below by a positive 
constant and Vuq G -L^(R^), smoother weak solutions may be built. Their existence stems from a 
quasi-conservation law involving the norm of Vu G L°°((0, T); L^(R^)) and of 9fn, Vp and V^n G 
L^((0, T); L-^(R^)) for any T < oo. For both types of weak solutions however, the problem of 
uniqueness has not been solved. 



Date: Dec. 22. 2012. 



1 



2 



M. PAICU, P. ZHANG, AND Z. ZHANG 



Ladyzenskaja and Solonnikov [16] first addressed the question of unique solvability of (1.1). More 
precisely, they considered the system (1.1) in a bounded domain $7 with homogeneous Dirichlet 

2—- V 

boundary condition for u. Under the assumption that uq & W p''^(Q) (p > d) is divergence free 
and vanishes on and that po G C^(0) is bounded away from zero, then they [16] proved 

• Global well-posedness in dimension d = 2; 

2— — 

• Local well-posedness in dimension d = 3. If in addition uq is small in W ^'^'(0), then 
global well-posedness holds true. 

More recently, Danchin [8] established the well-posedness of the system (1.1) in the whole space R*^ 
for small perturbations of some constant density. Abidi, Gui and Zhang [3] investigated the large 
time decay and global stability to any global smooth solutions of (1.1). 

Another important feature of (1.1) is the scaling invariant property: if {p,u) is a solution of 
(1.1) associated to the initial data (po)^o)) then [p{X'^t, Xx), Xu{X'^t, Xx)) is also a solution of (1.1) 
associated to the initial data {po{Xx) , Xuo{Xx)) . A functional space for the data {po,uo) or for the 
solution (p, u) is said to be at the scaling of the equation if its norm is invariant under the above 
transformation. In this framework, it has been stated in [1, 7] that for the initial data {po,uo) 
satisfying 

(po - 1) e bIi{k^), uo e bI~\k^) with divuo = 
and that for a small enough constant c 

\\P0 — 1|| d + \\uo\\ d_, < c, 

we have for any p G [1, 2d) 

. i 

• existence of global solution {p,u,Vp) with p — 1 G Cb([0, oo); ^^ ^^(R'')), u G C(,([0, oo); 

s|~\r'^)), and dtu,V^u,Vp G Li(R+; s|r'(R^)); 

• uniqueness in the above space if in addition p < d. 

These results have been somewhat extended in [2] so that uq belongs to a larger Besov space. Paicu 
and Zhang [18] further extended the well-posedness result in [2] so that even if one component of 
the initial velocity is large, (1.1) still has a unique global solution. The smallness assumption for the 
initial density in [1, 7] has also been removed in [4], and the restriction of p G [l,(i] for uniqueness 
result in [1, 7] has been removed recently in [9]. 

A byproduct in [4] implies the global existence of solutions to (1.1) in 3-D with initial density in 
L°°(R^) and having a positive lower bound, and initial velocity being sufficiently small in H^(It^). 
The authors [9] proved the global wellposedness of (1.1) provided that 

WPQ - 1|| d_l + ||«0ll d-l < C, 

d d -j^ 

for some sufficiently small constant c, where A4 {B^ ^ (R^) ) denotes the multiplier space of B^ ^ (R*^) . 
This space in particular includes initial densities having small jumps across a interface. 

Again in the scaling invariant framework, the authors [13] proved the global existence of weak 
solutions to (1.1) provided that the initial data satisfy the nonlinear smallness condition: 

(IIPo^ - + 11^0 II .-i+d) (^wiCrWuQf'' d) < CqH 

for some positive constants co,Cr and l<p<(i, l<r<oo, where Uq = (uq,--- ,Uq~^) and 
Uq = {uq,Uq). With a little bit more regularity assumption on the initial velocity, they [13] also 
proved the uniqueness of such solutions. 
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In general when po G L°°(R'^) with a positive lower bound and uq G H'^(R,'^), Danchin and 
Mucha [10] proved that the system (1.1) has a unique local solution. Furthermore, with the initial 
density fluctuation being sufficiently small, for any initial velocity uq G 5|2(H.^) H L^(R^) in two 

space dimensions, and uo G Bq^p'' (R'') with l<p<oo,(i<g<oo and 2 — | / |, they also proved 
the global well-posedness of (1.1). 

On the other hand, Hoff [14, 15] proved the global existence of small energy solutions to the 
isentropic compressible Navier-Stokes system. The main idea in [14, 15] is that with appropriate 
time weight (see Remark 1.2 for details), one can close the energy estimate for space derivatives of 
the velocity field even if the initial velocity only belongs to L^(R'^). Motivated by [14, 15], we shall 
investigate the global well-poscncss of (1.1) with less regular initial velocity than that in [10] and 
without the small fluctuation assumption on the initial density. We emphasize that the Lagrangian 
idea introduced in [9, 10] will also be essential for the proof of the uniqueness result here. 

Our main results in this paper can be listed as follows. 

Theorem 1.1. Let s > 0. Given the initial data {po,uo) satisfying 
(1.2) < Co < poix) < Co < +00, uo G H'{K^), 

the system (1.1) has a unique global solution (p, u) such that 

Co < p{t, x) < Cq for {t, x) G [0, +oo) x R^, 

Aoit)<C\\uo\\l2, 

Ai{t)<C\\uo\\hew{C\\uo\\l,}, 

A2{t) < C(l + \\uo\\l^){\\uo\\jjs + \\uo\\jjs)exp {CWuoWi^}, 

for any t G [0, +oo). Here C is a constant depending on co,Co, and Ao{t),Ai{t), and A2{t) axe 
defined by 

Ao(r)=^-sup [ p\u{t,x)fdx+ [ [ \Vufdxdt, 
2 te[o,T] J-R? Jo Jn.^ 

Ai(T)=^l- sup a{t)^-' [ \Vu(t,x)\'^dx+ C [ a{t)^-' {p\ut\'^ + \V'^u\'^ + \Vp\'^) dxdt, 

^ te[o,Ti Jr2 Jo Jr2 



A2(T)'i'l sup aitf-^ [ {p\ut{t,x)f + \VMt,x)f + \Vp\^)dx 

+ f f (T(t?~'Nut\'^ dxdt, 
Jo iR2 

with a{t) =^min(l,f). 

Theorem 1.2. Given the initial data {po,uo) satisfying 

(1.4) < Co < /9o(a;) < Co < +00, uqEH^R^), 
there exists a constant £o > depending only on Co such that if 

(1.5) ||'"o||l2||Vuo||l2 < £o, 

the system (1.1) has a unique global solution (p, u) which satisfies 
Co < p{t, x) < Co for {t, x) G [0, +oo) x R^, 

Bo{t) < WpIuoWI,, 



(1.6) 

D i4.\ ^ OllV7„. I, 

lL2) 

B^ij) < C(l + llVuolli^) llVnolli. exp{C(||uo|li2 + ||Vuo|li2)}, 



i?i(t) <2||Vno||2 
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for any t G [0, +00). Here C is a constant depending on Cq, and Bo{t),Bi{t), and B2{t) are defined 
by 

5o(r) =^ sup / p\u{t,x)f dx + 2 / \Vu\'^dxdt, 
te[o,T] Jr3 Jo jr3 

Bi(r)^^sup [ \Vu{t,x)fdx + 2 r f {p\ut\^ + \V^u\^ + \Vp\^)dxdt, 

t6[0,T] JR^ ^0 JR^ 

B2{T)=^ snv a{t) f {p\ut{t,x)\^ + \V'^u{t,x)\^ + \Vp\^) dx + I I a{t)\Vut\^dxdt. 
te[o,T] Jb? Jo Jn^ 

Remark 1.1. We should point out that we do not need the lower bound assumption for the initial 
density in the existence part of Theorem 1.2. Indeed, the constant C in (1.6) is independent of cq 
in (1.4). We can also prove the local existence and uniqueness solution to the system (1-1) even if 
the initial velocity does not satisfy the smallness condition (1.5). One may check Remark 2.1 for 
details. 

Remark 1.2. The powers to the weight a(t) in the energy functionals, Ai{t), Bi{t) for i = 1,2, 
are motivated by the following observation: let s be a negative real number and {p,r) G [l,oo]^, a 
constant C exists such that 

(1-7) C-'\\f\\^^ <||ll*-^e*^/IUHlL^(R+.ii)^^ll/llB^ ' 

(see Theorem, 2.34 of [6] for instance). In particular, if uq € -ff^^R^) for s G (0, l),Vuo G 
H'-^iR"^) ^ B^-^{B?) and V'^uq G B^-J^iR"^). Then according to (1.7), 

t^||e*'^Vuo||L2 < Clltxolli/^ and ||e*'^V^'Uo||L2 < CHuoIIj?". 

This in some sense explains the weights in (1.3). Similarly, we can explain the weights in (1.6). 

Remark 1.3. We should point out that we can not directly apply Theorem 1 of [10] concerning the 
uniqueness of solutions to (1.1) to conclude the uniqueness part of Theorem 1.1 and Theorem 1.2. 
Yet the Lagrangian idea in [9, 10] can be successfully applied to prove the uniqueness part of both 
Theorems 1.1 and 1.2. And the uniqueness result of Germain [11] can not be applied here either. 
The uniqueness result of [11] requires the density function satisfying Vp G L°°{[0,T]; L'''{R'^)), but 
here our density function only belongs to L°°{[0,T] x R*^). Moreover, the velocity field in Theorems 
1.1 and 1.2 does not satisfy the time growth condition in [11], especially in Theorem 1.1. 

2. Global solutions to (1.1) with large bounded density 

The purpose of this section is to present the proof to the existence part of both Theorem 1.1 
and Theorem 1.2. 

2.1. Existence of the solution in 2-D. 

Proof to the existence part of Theorem 1.1. Let be the standard Priedrich's moUifier. We define 

Po=je*Po, Uo=j€*UO- 

And we choose e so small that 

j<ph{x)<2Co, xeR\ 

With the initial data {pq, Uq), the system (1.1) in 2-D has a unique global smooth solution {p'^,u'^). 
In what follows, we shall only present uniform energy estimates (1.3) for the approximate solutions 
{p'^,u^). Then the existence part of Theorem 1.1 essentially follows from (1.3) for {p^,u^) and a 
standard compactness argument. The uniqueness part of Theorem 1.1 will be proved in Section 3. 
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To simplify the notations, we will omit the superscript e in what follows. First of all, applying 
the basic energy estimate to (1.1) gives 

(2.1) Mt) = l\\pluo\\l2<C\\uo\\l2 for teR+. 

While it follows from the transport equation of (1.1) and (1.2) that 

(2.2) Co < p{t, x) < Co for {t, x) G R+ x . 

To derive the estimate for Ai{t), we get, by taking the inner product of the momentum 
equation of (1.1) with ut, that 

/ plutl'^ dx + - [ |VuP(ix = — / p(u ■ Vu) • Ufdx, 
Jb? 2 dt 7r2 J^2 

from which, we infer 

Ai{t)=^ [ [ a{T)p\ut\'^dxdT + la{t) f \Vu{t,x)fdx 
^2 3) Jo Jr2 2 Jr2 

<- / / a{T)p{u-Vu) ■utdxdT+ / \Vu\'^ dxdr. 
Jo 7r2 Jo 7r2 

In this subsection, we shall frequently use the following version of Gagliardo-Nirenberg inequality: 

(2.4) II ^IIlp(r2) ^ ^ll'^lll2(R_2j||Va||^2|,^2j for 2 < p < oo. 
By virtue of (2.2) and (2.4), we obtain 

/ / a{T)p{u ■ Vu) • ut dxdr 
Jo Jr^ 

<C [ air)\\u{T)\\L4^u{T)\\L4p'^ut{T)h2dT 
Jo 

(2.5) <C [ a(T)||ix(r)|||2||V^x(r)|U2||A^.(r)|||2||p5„t(r)|U2dr 

Jo 

<Cs [ a(T)||u(r)||^2||V^x(r)||i.||Atz(T)|U2dr + 5 / a(r)||p^«t(T)||i. dr 
Jo Jo 

<CsAo{t) [\{T)\\Vu{T)\\l,dT + 5 f a{T)\\Au{T)\\l2dT + SMt), 
Jo Jo 

for any 6 > 0, where Cs is a positive constant so that C5 — )• oo as (5 — )■ 0. Whereas thanks to (1.1), 
we write 

(2.6) -Au + Vp = -p{ut + u-Vu), 

which along with the classical estimate on the Stokes system ensures that 

||V\||i2 + ||Vp|U2 <C7(||pnt||i2 + \\pu ■ Vu\\l2) 

1 1 

<C{\\pUt\\L2 + ||M||22||Vli||x,2||V^li|||^2), 

so that 

(2.7) I|v2n||i2 + ||Vp||l2 < C{\\p'2ut\\L2 + ||n||i2||V«||i2). 
Substituting (2.7) into (2.5) gives rise to 

(2.8) 1/ / a{T)p{u-Vu)-utdxdT\<CsAo{t) [ (T(r) ||Vu(r) 11^2 dr + C(5li(t). 
Jo Jb? Jo 
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Summing up (2.1), (2.3) and (2.7)-(2.8) and taking 6 sufficiently small, we obtain 

Mt)+ [ [ ^(TiT)\V\\'^dxdT <C(\\U0\\12 [ \\Vu{T)\\l2MT)dT+\\uo\\l2). 

Applying Gronwall's inequality gives 

(2.9) Mt)+ [ I cr(r)|V\|^dxdr < C||uo||i2exp{C||no||i2}. 

To obtain the estimate of ^i(t), we need to use an interpolation argument. For this, we consider 
the linear momentum equation 

p{dtv + u ■ Vv) — + Vp = 0, v(0, x) = vo{x). 

Then it follows from the same line to the proof of (2.9) and (2.7) that 

/ / {p\vt\^ + \'7^v\^ + \Vp\^)dxdT+ [ \Vv{t,x)\^dx<C\\vo\\jj,exp{C\\uo\\t2}, 
Jo JR^ Jn"^ 

[ [ a{T){p\vtf + \V^vf + \Vp\'^)dxdT + a{t) [ \Vvit,x)\'^ dx < C\\vo\\l2 exp {C\\uo\\l2}. 
Jo Jb? Jn^ 

We define the linear operator Tvq = Vv. The above inequalities tell us that 

||ri;o||L2 < Cexp{C||uo||^2}||'yo||i?i, 
||r?;o||L2 < Ca{t)'^ exp {C||uo|||,2}||vo||l2, 
from which and Riesz-Thorin interpolation theorem [12], we infer 

||Vt'||i2 = ||rt>o||L2 < C'cr(*)^^ exp {C||tto||^2}||uo||j?»- 

Further, we define a family of operators T^vq = a{t)i dtv{t,x) for Rez G [0, 1]. Then we have 

\\TiyVo\\L^{o,t;L^) < C exp {C||no ||^2 } 11^0 Hi/i , 
\\Ti+iyVo\\L^{o,t;L^) < exp {C||no ||^2 } || || L2 > 
for any y G R. Apply Stein interpolation theorem [12] to get 

||cr(T) V||p^^;j||^2||^2(o,t) = ll^l-s^^0||L2(0,t;L2) < exp {C||uo || ^2 } [jt'O • 

The other terms can be treated in a similar way. Therefore we have 

(2.10) Ai{t) < C\\uo\\h exp {C7||«o|lt2}. 

Finally, we manipulate the energy estimate for u. We first get, by taking the time derivative 
to the momentum equation of (1.1), that 

p{utt + u ■ Vut) - Aut + Vpt = -pt{ut + u ■ Vn) - put ■ Vn. 

Taking the inner product of the above equation with ut, and then using integration by parts, 
we write 



ld_ 
2di 



/ p\ut\^dx+ / jVutl^dx 

Jr2 7r2 

~ / PtWi]"^ dx — / pt{u ■ Vu) ■ utdx — / p{uf Vu) ■ Ufdx, 
Jb? Jb? vr2 
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from which, we infer 

= — a(T)^^^ / pAuA"^ dxdr — / (/(r)^^* / pAu ■ V u) ■ ut dxdr 

(2.11) ^0 7o Vr^ 

- / cr(T)^~* / p{ufVu)-UtdxdT + {2- s) / C7(r)^~"* / pl-Utpdardr 
Jo iR2 7o 7r2 

+ 5 + £ + D. 
It is obvious to check that 

(2.12) \D\<{2-s)A^{t), 
and 

ra(r)2-||p||Loc||Vn|U2||^.t|||4dT 
Jo 

(2.13) <C / a{Tf-'\\Vu\\LA\ut\\LA\^ut\\L2dT 

Jo 

<C a{Tf-'\\Vu\\l4p^ut\\l.dT + \A2{t). 
While noticing that pt = —u ■ Vp, we get by using integration by parts and (2.4) that 
A = -2 [ <T(r)2"* [ put ■ {u ■ Vut) dxdr 

Jo Jr2 

<2 / a{Tf-'\\u\\L4p\\L^\\Vut\\L4ML*dT 
Jo 

(2.14) <C I a(r)2"^||n|||,||V«|||,||ut|||,||Vut|||,dr 

Jo 

<C^*a(r)^-||n||i.||Vn||i.||«*||i.dr + ijJ|*a(r)2-^||Vt.t||i.dT 
<C||txo||i.^*<7(r)2-||Vt.||i.||tx*||i.dr + ii2(i). 
Along the same line, we write B as 

B = [ <T(r)^"* / p{u ■ Vu) ■ {u ■ Vut) dxdr 
Jo Jb? 

+ / airf'^ / pi(u ■ Vu) ■ Vn) • ut dxdr 
Jo Jr2 

+ / o-(r)^"* / p{{u (g) u) : V'^)u ■ ut dxdr 
Jo J-s? 



def 
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By virtue of Holder inequality and (2.4), one has 

Bi < [' a{T)^-'\\p\\L^\\u\\ls\\Vu\\L4Vuth2dT 
Jo 

<C I a(r)2-^||n|||,||Vn||i2||A«|||,||Vnt|U2dr 
<C airf-' \\uh2 \\Vu\\l, \\Auh2 dr + 
which along with (2.7) implies 

B, <c{j\{Tf-^\Vu\\Up'2Ut\\l,dT 

+ 1^ a{Tf-^ Ml. II Vulli. dr] + ^^A^{t) 
<C^*a(r)2-^||Vu||i.||;9^ut||i.dr + C||^o|li2Ai(i)2 + 1^2(0. 



The same argument gives rise to 
ft 
'0 



B2<j a(r)2-||p||i^||«|U4||V«||i8||p^Ut|U2dr 

<C [ a(r)2-^||^|||,||V«|U2||Au|||,||p5ut|U2dr 
Jo 

<C f a{Tf-muh2\\Vu\\l,\\p"2Ut\\l,dT+ fairf-^AuWl^dr 
Jo Jo 

<C [\iTf-'\\u\\L2\\Vu\\Up"2ut\\l,dT + C /V(r)2-(IIP^^tlli2 + ll^lli2||V«lli2 
Jo Jo 

<c(||no|U2 ^ a(r)2-^||Vt^||i2||p^ut||i2 dr + Ai{tf + ||uo||i2^i(t)') + ^Mt), 



and 

rt 



< I <T(r)2-||p|||^||«||ioo||V2n|U2||p^ut|U2dT 
Jo 

<C [ a{Tf-'\\u\\L4Au\\l4p'2uth2dT 
Jo 

<C [ a{Tf-mu\\L2\\p-2Ut\\hdT+ [ c7(r)2-^||n||i2||V«||i2||p^^xt|U2dr 
Jo Jo 

<c{\\uo\\hA,itf + ||no||i2^i(t)2 + 1^ airf-^VuWUp'^utWl dr) + ^^(t) 



'0 

Summing up the above estimates, we conclude that 



B<C{l + \\uoh2) fa{T)'-4Vu\\Uphut\\l,dT 
Jo 



(2.15) 

+ C{l + \\uo\\l.)A,{tf + ^A,{t). 
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Combining (2.11) with (2.12)-(2.15), we obtain 

Mt) <C{ (1 + WmWh) 1^ a{Tf-'\\Vu\\l2 Wp'^utWh dr 

+ {l + \\uo\\l.)Al{t) + Ai{t)]. 
Applying Gronwall's inequality and (2.10) leads to 

Mt) < C(l + ll^io|li2)(||no||l^. + ||uo|||f.)exp {CUt^oUla}, 
which together with (2.7) ensures that 

A2{t) <C{A2{t) + \\uo\\l,Alit)) 

<C(l + ||'Uo||i2)(||'Uo|||-» + ||'"o||ff<.)exp{C||uo||^2}, 
This together with (2.1) and (2.10) completes the proof of (1.3). □ 

2.2. Existence of the solution in 3-D. 

Proof to the existence part of Theorem 1.2. By mollifying the initial density po, wc deduce from [4] 
that (1.1) has a unique global solution {p^,u^) provided that eq is small enough in (1.5). Then the 
existence part of Theorem 1.2 follows from the uniform estimate (1.6) for {p^,u^) and a standard 
compactness argument. For simplicity, we only present the a priori estimates (1.6) for smooth 
enough solutions (p,u) of (1.1). The uniqueness of such solution will be proved in Section 3. As a 
convention in the rest of this section, we shall always denote by C a constant depending on Cq in 
(1.4), which may be different from line to line. 

First of all, it is easy to check from (1.1) and (1.4) that 

Co < p{t, x) < Co for (t, x) G [0, +oo) x R^, 

(2.16) r fT r r 

/ p\u{t,x)\'^ dx + 2 / / \Vu\'^ dxdt = / po\uo\'^ dx. 
jr3 jo jr^ 

While we get by taking the inner product of the momentum equation to (1.1) and ut that 

(2.17) 2 /" /" p\ut\'^dxdT+ f \Vu{t,x)fdx<\\Vuo\\l2-[ [ p{u ■ Vu) ■ Ufdxdr. 

Jo Jr3 7r3 Jo Jr3 

In what follows, we need to use Gagliardo-Nirenberg inequality 

(2-18) l|a|lLp(R3) <<^ll«llf2(R3)||Va||22(^3) for 2<p<6. 

By virtue of (2.18), one has 

/•* i 

p{u ■ V u) ■ Ut dxdr < I \\p\\loo\\u\\i(i\\'\/u\\]^i\\p'^ut\\L2 dr 

JR^ Jo 

VuWl^Wl^uWl^Wp^UtWL^ dr 



Jo 

<C \\Vu\\l2\\Au\\L2dT + ]- I \\p^ut\\l2dT. 

Jo 4 Jo 



Whereas it follows from the momentum equation of (1.1) and classical estimates on the Stokes 
system that 

+ I|Vp||l2 <C{\\p^ut\\L2 + \\u ■ Vu\\l2) 
(2.19) , I 

<C{\\p2ut\\L2 + \\Vu\\l2) + -\\V^u\\l2, 
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SO that 

f f p{u-Vu)-utdxdT<C [ \\Vu\\l2dr+'^ [ \\p^ut\\l2dT, 
Jo jr3 Jo ^ Jo 

from which, (2.17) and (2.19), we infer 

(2.20) 5i(t) < ||Vuo||i2 + C\\uo\\l2Bi{tf. 

Hence, as long as we choose Eq small enough in (1.6), we obtain the estimate for Bi{t) in (1.6). 
We now turn to the estimate of B2{t). Indeed along the same line to the proof of (2.11), we have 

Ja(t) / p\ut\'^dx+ f a{T) [ \Vut\^dxdT 

= — / (T(r) / pt\ut^ dxdr — / (t{t) / pt{u ■ V u) ■ ut dxdr 

(2.21) Jo Jo Jb:-" 

— / (t(t) / p{ut ■ Vu) ■ Ufdxdr + / / p\ut\'^ dxdr 
Jo Jr.^ Jo Jr^ 

E + F + G + H. 
It is obvious to observe that 

(2.22) H < Bi{t) <2\\Vuo\\l2. 

Whereas using pt = —u ■ Vp and integrating by parts, we get, by applying (2.18), that 

E = — 2 (T(r) / p{u ■ Vuf) ■ Ut dxdr 
Jo Jr3 

/■* i i 

<2 / a{T)\\u\\Loo\\p\\l^\\Vut\\L4p^iJ-t\\L^dT 

Jo 

ft 



(2.23) 



Notice that 



Jo 

<c£(||Vt.||i. + ||V2tx||i.)a(r)||p§u,||i.cir+i£a(r)||V«*||i.dr. 



G< [ a(r)||p|||^||V«|U3||t.i|U6||pWi|U2dr 
Jo 

<C f a{T)\\Vu\\l,\\VMh\Nn\\L4p'^ML^dr. 







(2.23) holds also for G. To deal with F in (2.21), we get, by using pt = —u • V p once again and 
integrating by parts, to write 

F = (7(r) / p{u ■ Vu) • {u ■ Vut) dxdr 
Jo Jr3 

+ / '^i'^) / P((^ ■ V'") ■ ■ ^* dxdr 
Jo Jr3 

+ / (^{r) I p{{u^u) ■.V^)u-utdxdT=^ F1 + F2 + F3. 
Jo Jr3 
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By virtue of (1.6) for Bi{t) and (2.18) , one has 

Fi < a{T)\\p\\Lo^\\u\\l6\\Vu\\L6\\Vut\\L2dT 
Jo 

<C / c7(r)||Vu|||2||V^u||i2||Vut||i2 dr 
Jo 

<cj^ a{r)\\Vu\\l4VMhdT + ^B2{t)<C\\Vm\\l2 + ^B2{t). 
Along the same line, we have 

F2< a{T)\\p\\loo\\u\\L6\\Vu\\le\\p2Ut\\L2dT 
Jo 

ft 1 

<C / (7(T)||Vn||i2 llV^nlllz ||p2uj||i2 c/r 
Jo 

<C f (||Vu||i2 + \\VMhHr){\\p'2ut\\l2 + W'^Mh) dT. 
Jo 

The same estimate holds for F3, as 

-^3 < / (^{r)\\p\\loo\\u\\loo\\V'^u\\L2\\p2Ut\\L2dT 

Jo 

ft , 

<C / a{T)\\Vu\\i^2\\V'^u\\l2\\p2Ut\\L2 dr. 
Jo 

Therefore, we conclude that 

rt 1 1 

(2.24) F<C / (II Vtx||i2 + ||V\||i2)a(r)(||/92ut|||2 + II Vanilla) dr + C||Vno||i2 + -52(t). 

Jo 4 

On the other hand, notice from (2.19) that 

a{t){\\vMmh + llVp|li2) < Ca{t){\\p'^utml2 + \\Vum%) < C{B2{t) + ||Vuo||i2), 
which together with (2.21)-(2.24) ensures that 

B2{t) < c{l\\\Vu\\l2 + ||V2n|||2)a(r)(||plnt||i2 + HV^nUi^) dr + ||V«o|li2 + ||V«o|li2}, 

applying Gronwall's inequality gives rise to the estimate of B2{t) in (1.6). This completes the proof 
of (1.6). □ 



Remark 2.1. Along the same line to the derivation of (2.20), we also get 

l_d 
2di^ 



|Vn(t)||^2 + / p\ut\'^ dx = — / p{u ■ V u) ■ lit dx 

<C\\Vu\\l2 + ]- I p\ut\^dx, 
2 Jb? 

which gives 

(2.25) \\Vu{t)\\l2 + / p\utfdx < C\\Vu\\ 

at jj^3 



|6 

Il2- 



Hence if the initial velocity uq does not satisfy (1.5), we deduce from (2.25) that there exists a 
positive time T so that 

l|VM|ii^(^2) + Ii«tiii|(i2) < c||vuoiii2. 
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With the above estimate, we can obtain the estimate Bi{t) and B2{t) for t <% as we did before. 
This implies the local existence of solutions to (1.1) in 3-D with large data. 



3. Uniqueness of the solution 

3.1. More regularity of the solutions. Before we present the proof to the uniqueness part of 
both Theorem 1.1 and Theorem 1.2, we need the following regularity results for the solutions of 
(1.1) obtained in Section 2. 

Lemma 3.1. Let (p, u, Vp) he the global solution of (1.1) obtained in Theorem 1.2. Then for any 
T G R^, one has 



L 



'^a(t)(||Au(t)||ie + ||Vp(t)||ie)dt<C, 





T 



r 1 

(3.1) / \\Vu{t)\\L^dt<Ciaax{T^,T2), 

Jo 

a(t)5||V«(i)||ioc dt<C, 



I 

Jo 



for some constant C depending only on co,Co in (1.4) and \\uo\\fji- 
Proof. We first get by taking div to (2.6) that 

Vp = V(-A)-^ dw{p{dtu + u ■ Vu)}, 
which along with (1.6) and (2.6) ensures that 

(3.2) \\Au{t)\\L^ + \\Vp{t)\\L^<C\\{dtu + u-Vu)(t)\\L'i for any gG(l,oo). 

However, it follows from Sobolev imbedding theorem that 

||(n-Vn)(t)||i6 < ||n(t)|U-||A«(i)||L2 <C||V«(t)|||,||An(t)|||„ 
from which, (1.6) and (3.2), we infer that 

^ a{t){\\Au{t)\\le + \\Vp\\le)dt 

<C{ r ait)\\Vut\\l.dt+ sup (||Vu(t)|U2CT(t)||An(t)||i2) T ||An(t)||i. dtj < C. 
^JO te[o,T] Jo 

This proves the first part of (3.1). Then we deduce from it, Gagliardo-Nirenberg inequality and 
(1.6) that 

I Au||22 ll^'^ll^e dt 



[ \\Vu{t)\\L^dt<C [ 

Jo Jo 

<||An||J,(^,)(^^a(t)||A«(i)||iadi)^^^a(t)-5dt)^ 



<c(^j a(t)-5dt)' < Cmax(T2,T5). 
Along the same line, we can also prove the estimate for Jq cr{t)2\\S/u{t)\\'j^^ dt. □ 
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The 2-D version of the above lemma is more comphcated, which we present as follows. 

Lemma 3.2. Let {p, u, Vp) be the global solution of (1.1) obtained in Theorem 1.1. Then for any 
T G R+ and a G [0, 1), one has 



I 



^a(i)i+"-^(||(9i«, A«)(i)f 2 + l|VKi)f 2 ) ^ ^' 

T 



law 



(3.3) [ ||Vu(t)||Locdt < Cmax(r5(TfeT,r5), 

/ (7(t)^-^||Vu(i)||iocdt < C, 
Jo 

where the constant C depends on a, co,Co in (1.2) and WuqWh"- 

Proof. Fist of all, for any /3 G (0, 1), wc deduce from Gagliardo-Nirenberg inequality that 

^ < C||a||}^2'^||Va||^2 for any aeH^iK^), 
which along with (1.3) ensures that for any a, ^, 7 G (0, 1) satisfying ^ + 7 = 1 + 0; 
||a(t)^(i+°)||(pn.Vn)(t)|| . 11^, 
<Co||a(t)^(^+")||«(i)||^^||V«(i)||^^|L| 

<Co sup (||n(t)||i;^a(t)^^||Vtx(i)||?OI|l|V^(t)|lir^(i)'^^l|VVi)lll2|L. 

tG[0,T] T 

< C(co,Co, ||?xo||//0- 

Along the same line, we obtain the same estimate for ||c7(t)^'^2 \\dtu\\ 2 IL2 ■ On the other hand, 

II ^ ' II £1-0: W-Uj, 

it follows from (3.2) that 

\\Au{t) II 2 + II Vp(i) II 2<\\{ut + u-Vu){t)\\ ., 

from which and (3.4), we obtain the first inequality of (3.3). 

Whereas we get, by using Gagliardo-Nirenberg inequality once again, that 

\\Vu{t)\\Lo^dt<C [ ||V?x(i)||2f^||A?/(i)||^ dt 

Jo 

<C\\Vu\\f^^,^\\ait)'^\\Aum^^^^^ 

which together with (1.3) ensures the second inequality of (3.3). Along the same line, we can prove 
the last inequality in (3.3). □ 

3.2. Lagrangian formulation. As in [9, 10], we shall apply Lagrangian approach to prove the 
the uniqueness of the solutions. We remark that even with (3.1) and (3.3), the solution of (1.1) 
obtained in Theorem 1.1 and Theorem 1.2 does not satisfy the assumptions required by Theorem 1 
of [10] concerning the uniqueness of solutions to (1.1). Fortunately, the idea used to prove Theorem 
1 of [10] can be successfully applied here. 

Let {p,u,p) be the solution of (1.1) obtained in Theorem 1.1 and Theorem 1.2. Then thanks to 
(3.1) and (3.3), we can define the trajectory X{t,y) of u{t,x) by 

dtX{t,y) = u{t,X{t,y)), X{0,y) = y, 
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which leads to the following relation between the Eulerian coordinates x and the Lagrangian coor- 
dinates y: 

(3.5) X{t, y)=y+ f u(r, X(r, y))dr. 

Jo 

Moreover, we deduce from (3.1) and (3.3) that we can take T small enough such that 

(3.6) r \\Vu{t)\\Loodt <\. 

Jo ^ 

Then for t <T, X{t, y) is invertible with respect to y variables, and we denote by Y{t, •) its inverse 

dcf 

mapping. Let v{t,y) = u{t,x) = u{t, X{t,y)). One has 

dtv{t. y) = dtuit, x) + u{t, x) ■ Vu(t, x), 
(3.7) 

dxiU^{t,x) = dy^v^{t,y)dxiy'' for x = X{t,y), y = Y{t,x). 
Let A{t,y) =^ {VX{t,y))-^ = VxY{t,x). So we have 

(3.8) Vxu{t,x) = A\t,x)Vyv{t,y) and div u{t, x) = div {A{t,y)v{t,y)). 
By the chain rule, we also have 

(3.9) dlYy {A-)= A^Vy. 

Here and in what follows, we always denote A^ the transpose matrix of A. 
As in [10], we denote 

V„=^^*-Vj;, div„ =^div(^-) and Au=^divuV„, 

ri{t,y)'^^%{t,X{t,y)), vit,y) "^^^ u{t, X{t,y)) and n(t, y) t%(t, X(t, y)). 

Notice that for any t > 0, the solution of (1.1) obtained in Theorem 1.1 and Theorem 1.2 satisfies 
the smoothness assumption of Proposition 2 in [10], so that (77,1;, VII) defined by (3.10) solves 



(3.11) 



f dtv = 0, 

rjdtv - AuV + V„n = 0, 
div„ V = 0, 
_ {v,v)\t=o = {po,uo), 



which is the Lagrangian formulation of (1.1). 

Now we transform the regularity information of the solution in the Eulerian coordinates into 
those in the Lagrangian coordinates. 

Lemma 3.3. Let {p,u,\/p) be the global solution of (1.1) obtained in Theorem 1.2 and (?7,t;, 11) 
be given by (3.10). Then for any t <T determined by (3.6), one has 

f rHmv,V\){T)\\l, + ||Vn(r)||i3) dr < C, 
Jo 

(3.12) ||V^||ic«(i3) + ^ \\VviT)\\LoodT<Ct-4, 



ft 

i„_ , .,,2 



/ TmVv{T)WUdr<C, 
Jo 



for some constant C depending only on co,Co in (1.4) and ||uo||i^i. 
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Proof. We first deduce from (3.5), (3.6) and (3.1) that 

(3.13) W^yXit, OIIl- < exp{ / \\V:cu{t)\\l^ dr] < 

Jo 

which together with (3.1) and (3.10) implies that 

/ ||Vv(t)||Loo dt < Ct-i and / A ||Vt;(r)|||cx. dr < C. 
Jo Jo 

Furthermore, thanks to det(^^^) = 1, and (1.6), (3.1), one has 

On the other hand, it follows from the proof of (3.13) that 

||V2x(i,-)||LP <exp{C f\\VMr)\\L^dT} T ||VVt, X(t, 0)1^^ 

Jo Jo 

<C ||V^ u{t, ■)\\lp dr for any p G [l,oo]. 
Jo 

In particular, if we take p = 3 in the above inequality and use (3.1) to get 



from which and (3.10), we infer 

<C{1 + \\t-^VIu\\l2^l^)) 
<C7(l + ||V^tx||i,(^,)||rW^«|l!2(^6))<^^- 
On the other hand, thanks to (3.6), we have for f < T 

(3.14) A(t,y) = V,y(t,x) = (/d + (V^X(i,y)-/d))-' = ^(-l)'^(/ Vj,n(t', X(i', y)) dt') , 

t=o 

for X = X{t, y), so that 

W^AU^m <c\\vlu{T,x{T,-mLim\\^yX\\Lr{L^) < ctl 

Finally, it follows from (1.6), (3.1) and (3.7) that 

<C[\\drU\\l,^^,^\\T-2VdrU\\l,^^^ 

This completes the proof of the lemma. 
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Lemma 3.4. Let {p,u,Vp) he the global solution of (1.1) obtained in Theorem 1.1 and (7^,^,11) 
he given by (3.10). Then for any t <T determined by (3.6) and < a < s, one has 

tr''"'-'{\\{drV,V\){T)f , + ||Vn(r)f . ) dr < C, 

ft 3 

(3.15) IIV^II 2 <Ct^, / \\Vv{t)\\loo dr < Ct^^ , 

ft 

-^~Tji:z II V7„,/'^M|2 



f T''^\\Vv{T)\\l^dT <C, 

Jo 



for some constant C depending only on co,Co in (1.2) and H'Wolli?"- 

Proof. The proof is similar to Lemma 3.3. We omit the details. □ 
3.3. The proof of the uniqueness. We first recall the following lemma from [10]. 

Lemma 3.5. Let rj G L°°(R'^) be a time independent positive function, and be bounded away 
from zero. Let R satisfy Rt G L'^{{0,T) x R'') andVdivR G L'^{{0,T) x R''). Then the following 
system 

' r]dtv -Av + Vn = f {t,x) € (0,T) x R'^, 
divv = div R, 
v\t=o = Vq, 

has a unique solution {v, VII) such that 

W^vh^iL^) + IK^t, V^^;, Vn)||^2 (^2) < C{\\Vvo\\l^ + \\{f,Rt)\\LUL^) + ||Vdivi?||^2,(^2)), 

where C depends on inf rj and supr/, but independent ofT. 

Proof to the uniqueness parts of Theorems 1.1 and 1.2. Let {pi,Ui,'Vpi),i = 1,2, be two solutions 
of (1.1) obtained in Theorem 1.1 and Theorem 1.2, and (r/j, Wj, VII), ? = 1,2, be determined by 

(3.10). We denote A A{ui) for i = 1, 2, and Sv =^ V2 - vi,6U =^ U2 - Hi and SA = A2 - Ai, 
then we deduce from (3.11) that 

PodtSv - A5v + VSU = Sfi + Sf2, 
(3.16) ^ div 5v = div dg, 

Sv\t=o = 0, 



with 



Sfi [(V - v„jni - (V - v„,)n2] 



(3.17) 



{Id- Al)V6U - 6A^VUi, 

6f2 [(A - A„Jt;i - (A - A,Jv2 

div[{Id - A2Al)V6v + {AiA\ - A2A*)Vi;i], 

- {Id - Ai)vi + {Id 
{Id- A2)Sv - SAvi. 



Sg =^ - {Id - Ai)vi + {Id - A2)v2 



We denote 



5E{t) -i^* ||V,5i;|Uoo(i2) + \\{dt5v,VHv,V5mq{L^y 
Then we infer from Lemma 3.5 and (3.16) that 

mt) < C^(I|5/i||l2(^2) + ||5/2||l?(l2) + ||Vdiv55||i2(^2) + ||5t55llL?(L2))- 
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We will show that 

(3.18) Wihf^L^) + ||^/2|L2(^.) + ||Vdiv^5llL?(L2) + ll^t^5llL?(L2) < e{t)5E{t), 

where the function £{t) tends to zero as t goes to zero. With (3.18) being granted, we infer that 

SE{t) < £{t)5E{t), 

which ensures the uniqueness of solutions obtained in Theorem 1.1 and Theorem 1.2 on a sufficiently 
small time interval [0, Ti]. The uniqueness on the whole time [0, oo) can be obtained by a bootstrap 
argument. □ 
Now let us turn to the proof (3.18). Indeed thanks to (3.6), we can take the time T to be small 
enough so that 

/ ||Vz;i(r)||L^dr < J, z = l,2. 
Jo ^ 

As a convention in the sequel, we shall always assume that t <T. Thanks to (3.14), we write 

(3.19) 5A{t) = {fv5vdT){Y,T.'^i'^2~^~') ^it^ Ci{t)^^^ f Vv^dr. 

The proof of (3.18) will split into the following two cases. 
Proof of (3.18) in 3-D case. We first deduce from (3.14) and (3.12) that 

||(/d-4)V5n||^2(i2) < C j\\VyV2{t')\\L^dt'\\V5Il\\L2^i^2) < CtUE{t). 

While it follows from (3.19) that 

WSAU^^i^e^ < C\\£\VSv\dt'\\^^^^,^ < Ct^\V^Svh2^^-2), 

which along with (3.12) implies 

||,5A*Vni 11^2(^2) <||r-i5^(T)|Uoo(i6)||rivni(r)||i2(^3) < CtUE{t), 
so that thanks to (3.17), we obtain 

(3.20) \\dfihuL^)<CtUE{t). 
Next we handle Vdiv^g^. We first get by applying (3.19) that 

(3.21) \V6A{t)\<C [ \V^6v\dT + C [ \VSv\dT [ (|VV| + |V\2|)dr. 

Jo Jo Jo 

Hence, we have by (3.12) that 
||V(MVi;i)|L2(^2) 

<C[\\£ \V^Sv\dr'Vvi\\^,^^,^ + \\J^ \VSv\dT'V\i\\^,^^,^ 
+ 11^" \V6v\ dr' + \V\2\)dT'Vvi\\^,^^,^'j 

+ \\r--* I |VHc^r'||^^^(^6)M(VV,V\2)|L2(^3)(l + ||TWi;i||^2(^^))} 
< CtUE{t). 
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Along the same line, one has 

iiv((/cz-^2)vmi^2(^2) 

< Cn5E{t). 
Thus, it follows from (3.9) and (3.17) that 

(3.22) \\VdiY5g\\^^L^^<Ct^5E{t). 
To deal with 5/2, we write 

(3.23) (^2^2 - AiA{)Vvi = {-5A{Id - A^) + {Id - Ai)dA^)Vvi + (M* + SA)Vvi. 

It is easy to check that 

||div[(M(/d-A*))Vt;i]||i2(i.) 

< ||r"3V(5yl||ico(i2)||/(i- ^2||L,°°(L«=)||r^Vvi||j;,2(icx=) 

+ lk~^'^^llL?°(L6)l|VA2||ioo(i3)||r2Vt;i||i2(ioo) 

+ ||r~5(5A||^oo(i6) \\Id - A2 llri V^vi ||l2(£,3) , 
which along with (3.21) and (3.12) implies that 
||div[(M(/d-4))Vi;i]|L2(^2) 

< Ci^ ||V2<5^;|L2(^2){||V7;2||L|(Loc)(||TiVi;i|U2(i^) + ||ri VV|L2(^3)) 
+ l|VA2||i^-(i3)||r2V?;i|U2(ioo)} < CtUE{t). 

The same estimate holds for div[{Id — Ai)dA^)Vvi\ and div[(M* + (5^)V'Ui] . Hence, 

II div[(A24 - AiA{)Vvi] 11^2(^2) < CtUE{t). 

To handle dw[{Id - A2A^2)'^^'^] > we write 

(3.24) {Id - A2Ai)VSv = -{{Id - A2){Id - A\) - {Id - A^) - {Id - A2)}V5v. 
Then we get, by using (3.12), that 

||div[(7d-^24)VHllLf(L2) 

< c{||V^2||L-(L3)||/d- ^2||Lr(L-)l|V5t;||i2(^6) + ||(/d- A2)||l,-(L-^ 

< Ct-i5E{t). 
Summing up, we obtain 

(3.25) 115/211^2(^2) <Ct35£;(t). 
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Finally, let us turn to the estimate of ||9t(5g'||^2(^2). Thanks to (3.19), we have 

\\dt[6Avi]\\^.^L-) <c{||(VM^i||l?(l2) + ||_^ \V5v\dT'\{VvuVv2)\vi\\^,^^^,^ 
+ \\j^^\V5v\dT'\dM\\^,^^,^] 
<c|||V5i;||ioo(i2)||i;i||i2(ioo) + l^'^^l'^'^'IL-(L6) 

X (||r3(VVl,VV2)||i2(^oo)||'Ul||^^^(^3^ + ||T5atVl||i2(i3))}, 

which together with (1.6) and (3.12) ensures that 

{1 1 J 

ll^^l|ll2(i2)||Aui 1122(^2) ||V5u||l^-(L2) +t4 II V25u||i2(^2) 

1 1 1 1 -| 

X (||t4(Vi;i, Vv2)||i2(i^)||ui||2c«(i2)||Vni||^^(^2) + \\TidtVi\\q^L^^)^ 

<Ct^\\V5v\\L^(^L2)+Ct^\\V'^5v\\L2^i^2) < CtUE{t). 
Along the same line, one has 

\\dt[{Id - ^2)HIIl?(L2) <c{\\Vv25vh2^L^) + 11^' \Vv2\ dr'\Vv2\\6v\\\^,^^,^ 

+ \VV2\dT'\dtSv\\\^,^^,^} 

<c{ll<^^^llL?°(L6)||Vi;2||Lf(L3) + l|Vi;2||Li(Loc) 

X (I|V^^2||l2(l3)II'^u||l-(L6) + \\dtSv\\L2(^L^))^ 

<CtUE{t). 

Hence, we arrive at 

(3.26) \\dt5gh2^L^2) < Ct\5E{t). 

Then (3.18) follows from (3.20), (3.22), (3.25) and (3.26). □ 

Proof of (3.18) in 2-D case. In what follows, we shall always take < a < s. By virtue of (3.14) 
and (3.15), we have 

||(/d- 4)VOT||i2(i2) < C^* \\VyV2{t')\\L^ dt'\\V6Jl\\Ll{L^) < CtW5E{t), 
and by Gagliardo-Nirenberg inequality, one has 

||M|| , 2, <C|| I \V5v\dt'\\ , 2, 

(3.27) <C \\V5v\\l2\\V^5v\\];^'^ dt' <Ct^6E{t), 

Jo 

which along with (3.15) implies that 

||M*Vni 11^2(^2) <||r-^M(r)||^^^^^|^||T^Vni(r)||^^^^^^ < CtUE{t). 
As a consequence, we obtain 

(3.28) ||5/i 11^2(^2) <Ct^5i?(t). 
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While due to (3.21), we get, by using (3.27) and (3.15), that 

\\V{SAVvi)h2^L.) < c{\\l^\V'Sv\dT'Vvi\\^,^^,^ + ||^V'^^MrW||^.(^,) 

+ 11^" \VSv\ dr' {\V\i\ + \V\2\)dT'Vvi\\^,^^,^] 

+ ||r-^^ [ \VSv\dT'\\^^^2 \\T^^{V^vi,V'^V2 

Jo y^") 

X (l + ||T^Vi;i|L2(^c.))} 

and we also have 

\\V{{Id-A2)VSv)h.^L,) 

<c{\\j\v\,\dr'l^^^^^\\VSv\\^,^^,^ 
Here we used the fact that 

(3.29) W^^^huLh - ^ CtfSE{t). 
Hence, we obtain 

(3.30) llVdiv 5511^2(^2) < CtW^<5£:(t). 

Now we handle ||(^/2||l2(l2)- Indeed it follows from (3.15) and (3.27) that 
||div[(,5^(/d-4))Vi;i]|L2(^2) 

< llr-ter V(5^||ioc(i2) \\Id - A2|Uoo(ioc) llrter Vt;i 11^2(^00) 

The same estimate holds for the remaining terms in(3.23). Hence, we get 

II dw[{A2Ai - AiA\)Vvi] 11^2(^2) < Ct^dE{t). 
Whereas thanks to (3.24), we get, by using (3.15), that 

iidiv((/d-A24)vmL2(^2) 

<C'(||VA2||^^^^^^^||/d-^2||L-(L-)||V5^;||^,^^|^ + ||(/d-^2)||Lc«(^^)||^ 

< Ct^(^6E{t). 
So we obtain 

(3.31) \\Sf2h2^i^2)<Ct^5Eit). 
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Finally we deal with \\dtSg\\j^2(^j^2y As in the 3-D case, we have 

\\dt[6Av,]\\L2^L^) <c{||(VH^^i||l|(l2) + \\£ \VSv\dr'\{Vvi,Vv2)\v,\\^,^^,^ 



+ 



\J^\VSv\dT'\dtV,\\^^^,^]. 



By (3.15), the first term on the right hand side is bounded by 

1 9 i s 

\NM\l^{l2)\\'"i\\lI(loo) < C6E{t)\\vi\\l^^^2^\\V^vi\\l,^^2^ < Ct^5E{t), 
and the third term is bounded by 

1+^ r 1^ 1^ ,|j II 1+1^^ II ^ <Cti6E(t), 

and the second term is bounded by 

l+g-s (l-a)a l+a-3 (l-s)ct 

< Ct 

On the other hand, it follows from Gagliardo-Nirenberg inequality that 
from which and (3.15), we infer that 

||5t[(/d-^2)HllL?(L2) 

< c{\\Vv2Svh2^L-^ + \\J^^\Vv2\dT'\Vv2\\6v\\\^,^^,^ + \\£\Vv2\dT'\dtSv\\\^,^^,^] 

<C<,\\t 2 Sv\\ , 2J|t 2 Vf2ll 2 

+ ||V^;2||Li(L^)(||T-^,5i;||^^^^^|^||r^V^;2||^,^^^^ + ll^f^HLf (l2)) } 

< Ct^5E{t). 
Therefore, we arrive at 

(3.32) \\dM\L2^L'2) < Ctf5E{t). 

Then (3.18) in the 2-D case follows from (3.28)-(3.32). □ 
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